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Abstract Existing theory and models suggest that a Type I (merger) GRB should have a 
larger jet beaming angle than a Type II (collapsar) GRB, but so far no statistical evidence 
is available to support this suggestion. In this paper, we obtain a sample of 37 beaming 
angles and calculate the probability that this is true. A correction is also devised to account 
for the scarcity of Type I GRBs in our sample. The probability is calculated to be 83% 
without the correction and 7 1 % with it. 
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1 INTRODUCTION 

There are two intrinsically different phenomena that give rise to gamma-ray bursts (GRBs). One is 
the merging of two compact objects, such as that of a neutron star and a black hole; the other is the 
core collapse of a massive star "collapsar"), such as the birth of a Type Ib/c Supernova. According 
to a classification scheme developed recently (Zhang, 2007 and Zhang et al, 2009), a GRB obtained 
from the former channel is defined as a Type I GRB, while one from the latter would be classified as a 
Type II GRB. This classification scheme, though more intrinsic than the classical short/hard vs long/soft 
categories, is not easy to carry out at the present stage, since many different criteria need to be applied 
in order to determine the progenitor of a GRB, few of which are decisive. One of these criteria is that 
Type I GRBs are usually short/hard ( Tgo < 2 s ), while Type II GRBs are usually long/soft ( Tgo > 2 s 
). This criteria, though supported by many individual cases, is not decisive, as any short GRB could 
theoretically have been a long GRB had it occurred at a high enough redshift. Of the criteria that are 
decisive, such as whether or not a GRB emanates gravitational radiation, most are impractical and the 
remainder is not applicable to a great majority of GRBs at the present stage. 

GRBs eject their energy in the form of jets. These jets have already been modeled, and their beaming 
angles (also called opening angles) can be calculated from observed physical data (Sari et al, 1999). The 
size of the beaming angle, according to the models, is highly dependent on the nature of the GRB 
progenitor: theoretically, Type II GRBs should have smaller beaming angles than Type I GRBs, due 
to the collimation effect of the stellar wind hugging collapsars. However, this dependency has not yet 
been supported statistically, possibly due to the aggravatingly small sample of Type I GRB jet beaming 
angles available, hence this paper, which specifically seeks to find statistical evidence of this theoretical 
relationship. 

This paper is divided into 5 sections, of which this introduction is the first. In the second, we will 
list the relevant observational data and use it to obtain the jet beaming angles of 37 GRBs. The third 
section will cover data reduction and statistical analysis. The results and discussion of the analysis will 
be presented in the fourth section, and the entire paper will be summarized in the fifth section, which is 
the conclusion. 
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2 COLLECTED DATA 

According to the literature (Sari et al, 1999), the jet beaming angle of a GRB could be obtained by 
Eq.(l): 

0. = Q57(-^-)i (±±±)-l( Ei "M )- i (ItL)i( n )i (i) 
J l lday j 1 2 ' l 10 53 erg j l 0.2 j l 0.1cm- 3j ' ( ' 

where 6j is the jet beaming angle mentioned above, tj is the time of the jet break, z is the redshift at 
which the GRB took place, Ei SO (-f) is the isotropic energy released by the burst, ?; 7 is the gamma-ray 
radiative efficiency of the jet and n is the density of the interstellar medium, tj, z and Ei SO (j) are 
observable parameters, shown in Table 1 (Appendix A). r; 7 and n, however, are not observable, but 
thankfully do not have a very significant effect on the value of 9j. Therefore, we set ?? 7 = 0.2 and 
n = 0.1 cm~ 3 . The jet beaming angles obtained from these data are presented in Appendix A and in 
Fig. 1 below. In Fig. 1, the bars filled in black, stripes and blue correspond to Sample I, Sample II and 
unclassified jet beaming angles (defined below) respectively. 




Fig. 1 Jet beaming angles: complete sample. Sample 
I is shown in black filling, while Sample II is shown 
with striped filling. Filled in blue are those GRBs that 
have not been positively identified as either Sample I or 
Sample II. The horizontal axis is Oj in degrees, while 
the vertical axis the number of opening angles that are 
of the specified magnitude. 



3 DATA REDUCTION 

As shown in Appendix A, we have only managed to confirm 2 Type I bursts and 1 3 Type II bursts among 
the GRBs whose parameters we obtained. 

Of the remaining GRBs, one is identified to be a short burst, while 21 of the others are known to be 
long. According to experience, short GRBs are usually of a Type I origin, while long bursts are usually 
of Type II. Therefore, we incorporate the short GRBs into the Type I sample, and the long GRBs into 
the Type II sample. It is hoped that, statistically, this will serve the purpose of giving a larger sample 
while maintaining sample integrity. 
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Six bursts, namely GRB970828, GRB990510, GRB990705, GRB991216, GRB000301C and 
GRBO 10222 cannot be identified as either Type I/1I or long/short bursts, and therefore cannot be used 
in the statistics. These bursts are represented by the columns filled in blue in Fig. 1. 

Thus, we obtain a "Type I + short" sample, noted as "Sample I" for the rest of the paper, and a 
"Type II + long" sample, hereby noted as "Sample II". Their jet beaming angles are presented in Fig. 
1 (black and striped filling respectively). A primary contributing factor to the scarcity of Sample I data 
is the fact that tj is hard to determine in short bursts, rendering Eq. (1) unapplicable. Of the five bursts 
known to be of Type I origin (Zhang et al, 2009), we obtained the opening angles for two of them (see 
Appendix A for details). Of the remaining three (GRB050509B, GRB050724 and GRB061006), the jet 
break time of GRB050509B and GRB06 1 006 cannot be found, while only a lower limit constraint could 
be placed on the jet break time of GRB050724. The problem of whether lower limit constraints could 
be used in the sample is discussed further in section 4. Of the short/hard bursts listed in Zhang et al, 
2009, the opening angle of only one could be ascertained. Since GRB080503, the latest short/hard burst 
mentioned in Zhang et al, 2009, five bursts (GRB080702A, GRB080905A, GRB080919, GRB081024 
and GRB090510) have been observed according to GCN reports that fall indisputably in the short/hard 
category. However, no jet break time has been found for any of them. 

From Fig. 1, it should be quite obvious that Sample I data have a greater arithmetic mean in com- 
parison to Sample II data. This is indeed the case, as Sample I has a mean value of 10.42 degrees, while 
the mean of Sample II is only 3.42 degrees. This is supportive of the statement that Type I GRBs have a 
larger beaming angle than Type II GRBs. However, the fact that there are 4 of the Sample II data that are 
larger than 2 of the 3 Sample II data leaves plenty of room for argument. Therefore, statistical analysis 
based on the data that lead to a quantitative result on the validity of the statement above is required. In 
order to obtain such a quantitative result, statistical fitting must be carried out. Many papers consider 
the Gaussian distribution under similar circumstances (for instance: Zhang and Meszaros, 2002), and 
therefore we will fit the data using the Gaussian distribution. Whether opening angles really do follow 
Gaussian distributions is indeed rather problematic: this topic will be discussed further below. 

Assuming that elements taken from Samples I and II are random variables that have Gaussian dis- 
tributions {N(fi, a 2 )), we take their probability densities to be 

A^) = ^=cxp[-i(^^) 2 ], (i = l,2), (2) 



respectively. Here, i = 1 corresponds to Sample I, while i = 2 corresponds to Sample II. fii and cr^ are 
the means and standard variations of the two samples respectively, and are hereby defined mathemati- 
cally as yUi = Xki and of = (Mi — Xki) 2 , where Xki is the "k"th datum from Sample i. 

k k 

Taking the sample mean and sample standard deviation as /i and a respectively for both samples 
(/ii = 10.42, H2 = 3.42, a\ = 46.24, o\ — 9.07), we obtain best-fit curves as shown in Figs. 2 and 
3, where they have been superimposed on their respective sample distributions for comparison. For the 
rest of this paper, these best-fit distributions for Sample I and Sample II will be termed Fit I and Fit II 
respectively. A Kolmogorov-Smirnov test has been preformed for Fit II, resulting in = 0.125, in 
other words, an acceptance level of 1 — a w 65%. Fit I is expected to have a very low acceptance level, 
which we did not calculate but are certain that it is (possibly much) smaller than 30%, which is one of 
the reasons why we decided to devise the correction below. 

Next, we calculate the probability that a random variable from Fit I (Oji) is larger than a random 
variable from Fit II (6* j2 )- The random variable 6j\ — 9j 2 should have a distribution of 

NiOji - 9 j2 ) = 1 exp[4( ^-^-^-^ f] ; (3) 

\Ja{ + cr|v27r ^ V c T + 2 

This distribution is another normal distribution with a mean of (/ii — /i 2 ) and a variance of (a 2 +cr|), 
as shown in Fig. 4. Thus, 
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Fig. 2 Sample I best fit probability distribution curve 
superimposed on Sample I bar graph. 
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Fig. 3 Sample II best fit probability distribution curve 
superimposed on Sample II bar graph. 



i.e. 



p(9 jl >e j2 ) = p(8 jl -e j2 >o)= / N(e 

Jo 



jl - V j2 )d{0jl - (7j2 



Pi 



> 



-(-oc 



1 



-_exp[--( 



1 X- (pi - yU 2 ) 



Calculating this by means of a Fortran program, we obtain 



) 2 }dx. 



(4) 



Pi{0ji > e j2 ) =0.83. 



(5) 
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Fig. 4 Distribution of random variable Oji — 9j2 - This 
is the probability distribution of the value of the differ- 
ence between a random variable from Fit I and another 
from Fit II. 

The probability shown in Eq. (5) is the probability that a Type I jet angle is larger than a Type II jet 
beaming angle, but only if Type I and II GRB beaming angles have a distribution strictly the same as 
Fit I and Fit II respectively. In the case of Fit II, 34 samples have been taken to make the fit, therefore 
the difference is considered negligible. However, for Fit I, which was made with only three samples, the 
effects of inconsistency must be considered. 

According to Eq. (4), P(9ji > 6^2) can be calculated once Hi and a\ are given as constants. 
But the hi and a\ in this equation have probability distributions of their own, which are reliant on the 
consistency of Fit I. They are 

w-v> (6) 

and 

respectively, where N — 3, the number of Sample I elements. Here, t and \ 2 correspond to a t distribu- 
tion and a \ 2 distribution respectively, X is the sample mean, Sx is the sample variance, Hi an d 01 are 
the theoretical mean and standard variation of Sample I respectively(shown here as variables). 

Integrating Eq. (4) over these probability distributions, taking ( Jy ~^) v ^ an( j ( N ~J) s x ^ y an( j 

z respectively, we write 

p 2 (e 3 i > M 

= J +oo X ^2)J^t(y 7 2)-Pidydz 
note that Pi, though given as a constant in Eq.(5), is dependent on (fii ~ M2) an d {p\ + ff). which are 
not constants for the purpose of the following calculations. Substituting the right hand side of Eq.(4) for 
Pi (y, z), and equivalent terms in y and z for \i\ an d we obtain 
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= J + °° X \ Zl 2) t(y, 2) • P 1 (y, z)dydz 




eMAC-M=^fW Z , (8) 



2 V 



where 2) = * (1 + X \z, 2) = \ex V {-\), f(y) = m = X - %y, and g(z) = a\ = 



(N-l)S 2 x /z. 

Again using a Fortran program, we calculate this final result to be 



P 2 (B jl >9 J2 ) = 0.71. (9) 



4 DISCUSSION 

Theoretically, it would have been more logical to perform a similar correction on Sample II in conjunc- 
tion with the one used on Sample I. However, this would lead to 5-dimensional integration, which would 
be far more than what our PC could manage in a short period of time (the 3-D integration for P2 took 
10 minutes at double precision, with a step length of 0.03 for x,y and z). Instead, we performed the 
3-D integration correction on Samples I and II respectively. Taking into account the first 4 digits after 
the decimal point, Pi is calculated to be 0.8272. This merely decreases to 0.8254 after the application 
of the correction to Sample II, while it decreases to P2 = 0.7118 after application of the correction to 
Sample I. From these numbers, we find that negligence of the Sample II correction induces an error of 
only the magnitude of 10~ 3 (i.e. the third digit after the decimal point), and thus find it safe to retain 
the first 2 digits without performing the correction to Sample II, hence the number of digits used in the 
presentation of our results (Eqs. (5) and (9)). 

With an acceptance level of only 65% for even Fit II, it is indeed debatable whether a Gaussian fit 
is appropriate for the Samples. However, due to the fact that there does not exist an established fit for 
opening angles, a Gaussian fit seems to be the natural choice. Further research may experiment on other 
parametrical fits that may yield better acceptance levels. 

It has been proposed that GRBs which would fall into the Sample I category, but have only minimum 
constraints for their jet break times (such as GRB050724) could be used in the statistics as well. In doing 
so, a minimum value could be calculated for both Pi and P2, on the basis of a larger sample. However, 
this could lead to complexities, since there are GRBs which would fall into the Sample II category 
that have only minimum constraints for their jet break times too. In this paper, we wish to avoid such 
complexities, and therefore use only exact data. 

The reader might want to notice that there are several problems that have not been taken into ac- 
count. Firstly, all data on Ei SO (^), z and tj in this paper have been collected from the literature. While 
this should cause no problems for z and Ei SO (7), it could be somewhat problematic for tj , since the time 
of the jet break is different for different bands at which the specific observation was made. Throughout 
this paper, this matter has been treated indiscriminately.Also, the data, as shown in Table 1 (Appendix 
A), are presented with a certain deviation in the literature. In other words, it is not precise. Even more 
prominent is the crisis that if we take Fit I to be exactly the distribution of Type I GRB jet beaming 
angles, a significant number of Type I jet beaming angles will be smaller than zero. These problems, 
throughout the paper, have been treated as insignificant details, hereby submitted to future scrutiny by 
the reader. 
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5 CONCLUSIONS 



In this paper, we have obtained a Sample I of 3 jet beaming angles and a Sample II of 34 jet beaming 
angles. These two samples are expected to be representative of Type I and Type II GRBs respectively. 
After that, normal (Gaussian) distributions have been fitted to the samples, resulting in Fit I and Fit EL 
Taking these fits to be representative of Type I and Type II GRBs, we then proceed to calculate the 
probability that a random variable taken from Fit I is larger than another taken from Fit H, thereby 
deriving the probability that a Type I GRB has a larger opening angle than a Type II GRB. Taking into 
account the uncertainty caused by the very small Sample I, we then devise a correction for the probability 
stated above, and reach a probability that is much lower but still significantly high nevertheless. 

If we take Sample I to be a sample perfectly representative of Type I GRBs, then it could be said, 
with an 83% degree of certainty, that Type I GRBs have larger jet beaming angles than Type II GRBs. 
However, this 83% drops to 71% once we take into account the uncertainty caused by our small sample 
in Sample I. In either case, it could be justifiably concluded, from our current samples, that Type I GRBs 
generally have a larger beaming angle in comparison to Type II GRBs. 

Our results are supportive of current models and theory. Type I GRBs in general have larger beaming 
angles in comparison to Type II GRBs, with a fairly high degree of certainty, though by far not high 
enough to become a decisive criterion as to whether a GRB is of Type I or Type II origin (so far). 

It is possible that with a larger Sample I, higher levels of certainty could be reached. Therefore, 
further observations that yield data concerning Sample I jet beaming angles (preferably Type I GRB 
jet beaming angles) are required for progress. It is also possible that with a large enough Sample I, the 
correction methods used in this paper (i.e. the correction for Sample I) will become obsolete, however, 
given the current rate at which Sample I opening angles are being derived, it seems unlikely that this 
would be the case in the near future. Lastly, it would be desirable to find a parametrical fit which could 
be supported theoretically, and which yields a higher level of acceptance than the Gaussian. 
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Appendix A: OBSERVED DATA & DERIVED JET BEAMING ANGLES 

Table 1 Data from Observations & Derived Data 
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GRB 




t 


z 


Category 


Oj(deg) 


050709 


0.00069' 1 ! 


lO^ 


0.161 1 ) 


Type 


18.1957 


060614 


0.02ll 2 l 


1.3W 


0.13W 


Type lI 2 J 


5.57885 


970508 


0.06ll 2 l 


25W 


0.835W 


Type III 2 ! 


12.3366 


980703 


0.72 M 


2.49I 4 ! 


0.966W 


Type III 2 ! 


3.71802 


990123 


22.9^ 


1.8W 


1.6 W 


Type III 2 ! 


1.92367 


990712 


0.0672I 5 ! 


1.6I 5 ! 


0.433I 5 ! 


Type III 2 ! 


4.77012 


000418 


0.75137I 6 ! 




0.1181I 6 ] 


Type III 2 ! 


8.54109 


000926 


2.7ll 2 l 


2.03W 


2.307W 


Type III 2 ! 


2.40104 


011211 


0.67234I 6 ! 


1.77(6] 


2.1416] 


Type IlPl 


2.76815 


020405 


IP] 


2. 74 W 


0.689 [41 


Type III 2 ! 


3.91555 


020813 


6.6® 


0.46W 


1.254W 


Type III 2 ! 


1.42145 


000328 


3.61I 5 ! 


0.8l 5 l 


1.52I 5 ! 


Type III 2 ! 


1.80902 


030329 


0.166I 5 ! 


0.5I 5 ! 


0.169I 5 ! 


Type III 2 ! 


2.97281 


041006 


0.83 [5 1 


0.16 [5] 


0.716 [5] 


Type nPl 


1.37314 


050525A 


0.25P1 


0.16W 


0.606W 


Type III 2 ! 


1.63548 


051221A 


0.024I 7 ! 


4.1W 


0.5465W 


short! 8 ] 


7.50341 


010921 


0.13611I 6 ] 


33^ 


0.4509I 6 ! 


long! 9 ! 


13.5235 


020124 


2.15^ 


3 [B] 


3.198I 5 ! 


long! 9 ! 


2.61651 


021004 


0.55601I 6 ! 


7.6l 6 l 


2.332I 6 ! 


long PI 


4.78799 


030226 


0.67^ 


0.84I 5 ! 


1.986I 5 ! 


long! 9 ! 


2.13395 


030429 


0.173I 5 ! 


1.77I 5 ! 


2.656I 5 ! 


long! 9 ! 


3.09816 


050315 


0.49I 10 ! 


2.6W 


1.95W 


long! 11 ! 


3.40529 


050318 


0.22I 3 ! 


0.12W 


1.44W 


long! 11 ! 


1.27531 


050401 


5.323I 12 ! 


0.06I 12 ! 


2 . 9 [12] 


long! 11 ! 


0.55384 


0504 16A 


0.0083I 5 ! 


its] 


0.653I 5 ! 


long! 11 ! 


4.92335 


050820A 


9.74(3] 


3.99W 


2.6lW 


long I 11 ! 


2.55110 


060124 


4.1 PI 


0.61W 


2.3W 


long I 11 ! 


1.45364 


060206 


0.43 [31 


0.82W 


4.05W 


long I 11 ! 


1.83552 


060210 


4.1 5 [13] 


2.16W 


3.91W 


long! 11 ! 


2.00912 


060526 


0.26I 3 ] 


0.98W 


3.2lW 


long! 11 ! 


2.23733 


060605 


0.25l 14 l 


0.27 [14] 


3.773 [141 


longf 11 ' 


1.32273 


060814 


0.7I 3 ! 


0.79W 


0.84W 


long I 11 ! 


2.48693 


061121 


2.25 [31 


0.28W 


1.31W 


long I 11 ! 


1.33753 


061126 


1.06I 15 ! 


25.7I 15 ! 


1.1 588 1 15 J 


long! 11 ! 


8.20785 


070125 


10.6 [161 


3.8l 16 l 


1.547I 16 ! 


long! 16 ! 


2.82483 


071010A 


0.036I 17 ! 


l[17] 


0.98l 17 l 


long! 11 ! 


3.83010 


0803 19B 


1 3 [18] 


0.032I 18 ] 


0.937I 18 ! 


long! 11 ! 


0.50876 


970828 


2.1982^ 


2.2I 6 ! 


0.9578t 6 ] 




3.09196 


990510 


1. 76349 ^ 


1.2^ 


1.6187I 6 ! 




2.27045 


990705 


2.55952I 6 ! 


1(6] 


0.8424t 6 ] 




2.30918 


991216 


5.35369I 6 ! 


1. 2 [6] 


1.02I 6 ! 




2.17822 


00030 1C 


0.43749I 61 


7.3l 6 l 


2.0335 [6] 




5.03377 


010222 


8.57841I 6 ! 


0.93I 6 ! 


1.4768I 6 ! 




1.72899 



References: [1] Fox et al. (2005); [2] Zhang et al. (2009); [3] Amati et al. (2008); [4] Liang et al. (2008); 
[5] Ghirlanda et al. (2007); [6] Bloom et al. (2008); [7] Soderberg et al. (2006); [8] Sakamoto et al. (2007); 
[9] Pelangeon et al. (2008); [10] Amati L. (2007); [11] Evans et al. (2009); [12] Kamble et al. (2009); [13] 
Ghirlanda et al. (2008); [14] Ferrero et al. (2009); [15] Perley et al. (2007); [16] Chandra et al. (2008); [17] 
Covino et al. (2008); [18] Racusin et al. (2008). 
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